Simulation of seismic response in a city-like 

environment 



Jean-Philippe Groby* Chrysoula Tsogka^ Armand Wirgin^ 

February 2, 2008 



Abstract 

We study the seismic response of idealized 2D cities, constituted by non equally-spaced, 
non equally-sized homogenized blocks anchored in a soft layer overlying a hard half space. 
The blocks and soft layer are occupied by dissipative media. To simulate such response, 
we use an approximation of the viscoelastic modulus by a low-order rational function of 
frequency and incorporate this approximation into a first-order-in-time scheme. Our re- 
sults display spatially-variable, strong, long-duration responses inside the blocks and on the 
ground, which qualitatively match the responses observed in some earthquake-prone cities 
of Mexico, France, the USA, etc. 
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1 Introduction 

A noticeable feature of many earthquake-prone cities such as Mexico- City [2], P], Nice 
0, Los Angeles ^ etc., is that they are partially or wholly built on soft soil. Seismogram 
records in such cities indicate amplification of the ground motion, beating phenomena, long 
codas and substantial spatial variability of response. 

To analyze the possible causes of these puzzling effects, we study the action of a seismic 

wave on a relatively-simple structural model with both geological and man-made features. Our 
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2D model has three components, as in [Hj, (from bottom to top in Fig. a hard half space 
(HHS), overlain by a soft dissipative soil layer (SL), in which are partially imbedded a set of 
even softer dissipative blocks (SB). HHS and SL are geological features, the set of SB, which are 
homogenized for purpose of the analysis, is man-made and constitutes the visible component of 
an idealized city. 

Real media disperse and attenuate waves. In [H], all the material components of the 'city' 
were taken to be lossless. In order to obtain a more realistic picture of wave propagation in 
the soft layer and in the homogenized blocks, we consider them to be composed of dissipative 
media. We follow the approach exposed in [HI; IHl and show how it can be used to simulate 
SH viscoelastic wave propagation. 

Even when the media are dissipative, we find that the presence of the blocks gives rise 
to ground motion amplification, beating phenomena, long duration and spatial variability of 
response, which are not present in the case of a flat ground (i.e., without homogenized blocks) 
for the same solicitation. We also find that these effects are similar, although with less amplitude, 
to those for non-dissipative media (as assumed in 0). 

In order to understand, identify, and quantify some causes of these features, we study the 
influence of the distance separating the homogenized blocks. If the homogenized blocks are 
considered to be homogenized buildings, we show that decreasing the distance between buildings 
increases the interaction between buildings, and between buildings and the soft layer, but does 
not modify fundamentally the response in terms of amplification and duration. We attribute 
duration effects to the excitation of Love-like modes due essentially to the presence of the blocks 

2 Basic Ingredients of our Model 

Referring once again to Fig. El the blocks of the city are assumed to be in welded contact with 
the substratum by an interface on which we assume (as in jO]) the continuity of displacement 
and normal stress. Our "city" is invariant in the Xs-direction with xi, X2, X3 being the cartesian 
coordinates, and X2 increasing with the depth (see FigEl wherein the sagittal plane is displayed). 
The support of the seismic source is a line in the xa-direction, and is located deep in the 
HHS, radiating a Ricker pulse cylindrical shear-horizontal (SH) displacement field. Only the 
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j;3-component of this field is non- vanishing and invariant with respect to X3, so that the total 
field underneath and on the free surface is also SH-polarized and invariant with respect to x^. 
The resulting problem is 2D with the displacement field depending only on xi, X2 and on time 
t. 

We denote by hi, Wi and j+i, the height, width of the block Bi and space interval between 
the blocks Bi and Sj+i for the three configurations, ; j = 1, 2, 3, we studied. The substratum 

is occupied by a linear, isotropic, medium Mq, characterized by mass density p{x) = po, shear 
modulus /i(x) = fiQ, and quality factor (5(x) = Qq = +00, with x = (xi, X2). The soft layer ili is 
occupied by a linear, isotropic, dissipative medium Mi characterized by mass density p(x) = pi, 
relaxed shear modulus Pr{x) = pm, and quality factor Q(x) = Qi. The blocks O2 are occupied 
by the linear, isotropic, dissipative medium M2 characterized by mass density /9(x) = p2, relaxed 
shear modulus Pr{x) = //jj2, and quality factor (5(x) = Q2. 

3 Methods 

Since problems in seismology are essentially concerned with transient phenomena it is natural 
to treat them in the space-time rather than the space-frequency framework. This also presents 
certain advantages from the numerical point of view. The obvious approach is then to solve 
the second-order-in-time wave equation (Navier equation when the medium is isotropic) for the 
displacement vector. A less-obvious, but numerically-advantageous method [l2], [IH], [H], is to 
cast the wave equation into the form of two first-order-in-time partial differential equations, for 
the two unknowns constituted by the displacement and pseudo- velocity vectors. 

In the following, we show briefiy (more details are given in 

0, 0, ini) how to do this 

for motion in viscoelastic, isotropic media. Since our problem relates to 2D, SH motion, the 
unknowns reduce to the single non-zero component of the displacement vector and the two 
components of a pseudo- velocity vector. We show that these unknowns satisfy two coupled first- 
order-in-time partial differential equations, in which appear a set of unknowns (usually called 
memory variables) whose existence is due to the viscoelastic nature of the media. Each of these 
coefficients is shown to satisfy a first-order-in-time partial differential equation whose driving 
term is related to the pseudo-velocity vector. 

In an isotropic, non pre-stressed, elastic ( viscoelastic) medium M, the stress tensor s^i is 
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related to the strain tensor e^i by the Hooke-Cauchy constitutive relation 

Ski = >^Skiemm + '^'fJ-eki ; k,l = 1,2,3. (1) 

wherein the Einstein summation convention is implicit, 5ki is the Kronecker delta and Cki is 
related to the particle displacement vector (whose elements are Uk) by 

^kl = ^ iuk,i + Ul^k) , (2) 

with Ufe,/ signifying the partial derivative of with respect to the Z-th cartesian coordinate xi. 
Furthermore, A is the bulk modulus and /x the rigidity. In elastic media, A and /x are, at most, 
functions of position, whereas in viscoelastic media (i.e., with memory), ^ depends also on time. 

The tensorial wave equation in elastic and viscoelastic media, in the context of linear (visco) 
elasticity, is 

Ski,k - pdfUi = -pfi ; k,l = 1,2,3 , (3) 

wherein p is the mass density of the medium, /; the l-th component of the applied force density 
vector, t the time variable and 5| := This set of partial differential equations can be 

recognized to be linear in terms of the displacement ui since Ski is a linear function of u;. 

The 2D SH aspect of the problem is encompassed in the assumptions: /i = /2 = 0, ui = 
U2 = and ^3^3 = 0. so that the tensorial wave equation becomes 

(Sl3,l + S23,2) - pdfUs = -pfs . (4) 

Let 

(7 := (si3, S23) , U--U3 , f := h . (5) 

Then 

(T = (AtW3,l, A*W3,2) = AiVu , (6) 

which is the expression of the constitutive relation for the situation corresponding to that of 2D 
SH motion in (visco) elastic media. 

By the same token, the wave equation for 2D SH motion in isotropic, non pre-stressed, 
(visco) elastic media becomes 

V-{pVu)-pd'iu = -pf , (7) 

which is of the general form 

V • {aVU) - hd^U + = . (8) 
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We are able to show that this equation can be obtained from the two first-order-in-time 
equations 

-dtV + aVU = , (9) 



hdtU -V -Y = cG . 



(10) 



by taking the divergence of the partial time derivative of (fTn|l . and subtracting the equations 
so obtained. This leads to provided G is chosen so that 



Thus, instead of being confronted with solving a single, second-order-in-time partial differen- 
tial equation for [/, we must solve two first-order-in-time partial differential equations for the 
displacement U and the pseudo-velocity V. 

By a viscoelastic solid we mean a material whose rigidity at an instant t depends on the 
whole strain history of the solid, so that the Hooke-Cauchy relation takes the more general 
(convolution) form (for 2D, SH motion) 



wherein x := {xi^X2)- When the constitutive relation ifT^ is inserted into the wave equation, the 
latter becomes an integro-differential equation which is very cumbersome to solve numerically 
because it requires saving in memory the whole history of the solution at all points of the 
compuational domain. To overcome this inconvenience, Emmerich and Korn |15| propose a 
(essentially Maxwell body) method whereby they approximate the Fourier spectrum /i(x, w) 
(with UJ the angular frequency) of the viscoelastic rigidity modulus /u(x, t) by a rational function 
of frequency. 

Consider .F(x, t) to be a function of x and t; then ^(x, w) is its Fourier spectrum (with oj 
the angular frequency) such that 



F = dtG . 



(11) 




(12) 



.F(x, t)= / J'^(x, w) exp(— iu;t)(i( 



(13) 



J —oo 



Consequently, the frequency domain constitutive relation corresponding to lfT2)l is 



cr(x, w) = /i(x, a;)Vii(x, w) . 



(14) 
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We suppose that the spectrum function fi{x, u) is known (from experimental data) and wish to 
approximate it in the rational function manner. To do this, we introduce the relaxation function 
i?(x, t) defined by 



Jt 



(15) 



This relaxation function is assumed in ^H] to take the form of a discrete sum of sinusoidal 
functions 



R{x,t) f» /Uk(x) 



Hit) 



(16) 



wherein J is some integer, preferably not too large to reduce the computational effort, H{t) is the 
Heaviside function, /i_R(x) the relaxed rigidity modulus (equal to the time-invariant rigidity in the 
elastic case), 6fi{x) a differential rigidity which is a measure of the departure from elasticity, ujj 
the relaxation frequencies, it being understood that the coefficients aj obey the (normalization) 
relation 

J 

" 1 . (17) 



J 

^aj(x) 



From JTSl-linilwe find 



/i(x, 



wherein 



^^(x) 



(19) 



In the works of Emmerich and Korn and Groby and Tsogka [Ij , [HI , (HI the weight functions 
{yj(x ; j = 1,2, .., J)} are determined in unique manner from a given ^{x,uj) by solving an 
overdetermined linear system arising from sampling (fT8|l at a set of relaxation frequencies 
which are chosen equidistant on a logarithmic scale in the interval [u; max /WO, u; max], with Umax 
the maximal frequency of the solicitation spectrum. 
We introduce the new variables ^^(x, w): 



Cj(x,u;) : = 



lUJ — IjJ 



■2/i(x)Mij(x)Vu(x, 



so that 



J J 
(T{x,t) = HRVu{x,t) +^Cj(x,t) = crij(x,t) +^Cj(x,t) • 



(20) 



(21) 
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It is then straightforward to show that Cj satisfies the differential equation: 

dtCji^, t) + cOjCji^, t) = yj{^)dt(TR{j^, t) . (22) 
We now introduce a new function r]j{-x.,t) such that 

p(x)a,\(x,t) := V-C,(x,t) . (23) 

Then, it can be shown that: 

d^Vji^, t) + ujjdtT]j{^, t) = -^V • (y,(x)/ii?(x)Vn(x, t)) , (24) 

J 

V • (^ij(x) Vn(x, t)) + /5(x) J2 9tVji^, t) - p(x)a2u(x, t) + p(x)/(x, t) = . (25) 

i=i 

The latter equation is of the form (jHl provided we make the associations: U = u, a = ^r, b = p, 
F = f + J2j=i^'t'nji ^iid c = /9, so that employing the previously-demonstrated equivalence 
between ^ and a system of two first-order-in-time partial differential equations, we find, at 
present, that is equivalent to the couple of equations: 

/iij(x)V^/(x, t) - atV(x, t) = 0, (26) 



J 

pdtu{^, t) - V • V(x, t) = pg{^, t)+Y, dtrjj (x, t) , (27) 

i=i 

wherein g is such that 

dtg{^,t) = f{^,t) . (28) 

The last step is to introduce into followed by integration over t and neglect of the 
integration constant: 

dtvji^, t) + ^jVji^, t) = • (y,(x)V(x, t)) . (29) 

Our procedure for solving a problem of 2D SH wave motion, in response to the solicitation 
/(x, t), in the temporal interval [0, T] and spatial domain Q, occupied by a isotropic, viscoelastic 
medium M characterized by the density p(x) and rigidity spectrum function p{x,uj), thus boils 
down to: 
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obtain g{x,t) from /(x, t) via 



dtgi^,t) = f{^,t) ; xgJ] , t € [0,T] , (30) 
obtain the weight functions yj(x) from Q := ^{fi{x,uj))/Q{n{x,uj)) via 



/i(x,a;) = /ii?(x) 



J 



j = l 3 



X G , u; G 



(31) 



solve the coupled system of three first-order-in-time partial differential equations 

fiR{x)Vu{x,t) - dtY{x,t) = ; x G , t G [0, T] , (32) 



pdM^, t) - V • V(x, t) = pg{x, t)+Y^ dtVji^, t) ; X G , t G [0, T] , (33) 

5t7?,(x, t) + iOjrjji^, t) = -^V • (?/j(x)V(x, t)) ; x G Q , t G [0, T] , (34) 
for the three unknown functions u{^,t)^ V(x, t), and ?7 := {?7j(x, t) ; j = 1, 2, .., J}. 



3.1 Description of the cities and methods of analysis 

The densities in the bedrock, soft layer and blocks (together with their foundations) where 
chosen to be: 2000fcg.m^^, 1300%. m^^, and 325kg. respectively. The bulk shear wave 
velocities in these three media were taken to be 600m. 60m. and 100m. s^^ (instead 
of 200m. s^^ in [6j) respectively, and the quality factors were chosen to be +oo, 30, and 100 
respectively. The foundation depth of the blocks was 10m and the soft layer thickness was taken 
to be 50m. The block widths and heights ranged over 30 — 60m and 50 — 70m respectively (see 
Table The block separations ranged over 60 - 100m for C\ 30 - 50m for and 10 - 30m 
for (see Table [TJ. 

We computed the seismic response both in the absence and presence of the blocks. The media 
below ground level were the same in both of these cases. Most of the aforementioned parameters 
are close to those of fTHI and and are fairly representative of the blocks and substratum 
at downtown sites in Mexico City. The computational domain was a 3500m x 3500m square 
discretized by a grid of 700 nodes in each dimension. This domain was surrounded by a PML 
layer 29 nodes thick, and 943 nodes where placed on the free surface |12| . 



To give a measure of the vulnerability of the blocks of the city, we computed the vulnerability 
index Rj, introduced in [6j: 



= rT . ' I2 (35) 



/o \dtu{yLj,t)\' 
/o^|a,nO(xo,t)r 

wherein: T is the time interval of significant shaking (taken herein to be 200s), \dtu{xj,t)\ the 
modulus-squared particle velocity at the midpoint of the j-th block, and \dtu^{xo,t)\'^ the same 
quantity measured on the ground in the absence of all the blocks. 

To give a measure of the strength of ground shaking, we computed the index i^jj+i: 

\dtu{xj^j+i,t)\^ 

Jo \dtu^ixo,t)\ 

wherein: T is as previously, \dtu{xjj-^-i,t)\^ is the modulus-squared particle velocity at the center 
of the ground segment between the the j — th block and the j + 1 — th block, and \dtvP{xQ,t)^ 
is as previously. 

3.2 Validation of the method 

We validated the time domain method on the canonical example of flat ground with no 
blocks. The viscoelastic modulus of the soft layer is given by Kjartansson's formula |17| : 



■ arctan^Qj^ ^ ) 

/"i(x,t^) = |/ii?(x,u;/j)| { — ^ ) (37) 



2 

-lui ' 



The source position was X5 = (Om, 3000m), and we computed the displacement at the point just 
below the hypothetic building Bi. In FigEl the small differences between the semi-analytical 
[TH] and numerical solution are due to: i) the fact that the semi-analytical solution involves 
numerical integration, and ii) discretization errors in our time-domain method. 



4 Results 

The aim of this work was first to prove that amplification, beatings, long duration, and 
spatial variability of response are observed even when the media (in the layer and blocks) are 
dissipative, and that essentially the same phenomena are involved in both the dissipative and 
non-dissipative cases. 

The snapshot of the modulus of the total displacement field in Fig. [HI pertains to the case 
in which there are no blocks in the city, for the dissipative layer case. As mentioned in ^Hl; [Eji 
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|18| . one does not expect the Love modes to be excited to any great extent when the source is 
far from the layer. This is what is actually observed in Fig. [HI since the displacement field in the 
soft layer is not typical of that of a Love mode. 

Fig. ini depicts the snapshot for a city overlying a dissipative layer with ten blocks. One 
observes a series of hot spots inside the low- velocity layer, constituting an indication of something 
like a standing wave in the layer betraying the excitation of a quasi-Love mode jTJ, 0], jTHj. 
This configuration also gives rise to rather large response in the blocks and on the ground, as 
was previously observed in [6j for non-dissipative media, and in for a periodic distribution 
of identical blocks. 

The time records and spectra, presented in Figs. [flHl call for the following comments. 

• Fig. El exhibits the time history of the velocity at the summit of building no. 6, for both 
the dissipative and non-dissipative cases. The durations are much longer when there is no 
dissipation. Nevertheless, in the dissipative case, beatings are observed and the duration 
of motion is of the order of 2 min. 

• Fig. ini contains the time histories of the velocity, and the velocity spectra at the summits 
and on the ground in between blocks for city C^. 

• Fig- El contains the time histories of the velocity, and the velocity spectra at the summits 
and on the ground in between blocks for city C^. 

• Fig. IHl contains the time histories of the velocity, and the velocity spectra at the summits 
and on the ground in between blocks for city C^. 

The results in terms of the vulnerability indices are given in Table 121 

5 Discussion 

This work was initiated in [HI, in which the authors considered the soft layer and blocks to 
be occupied by non-dissipative media. 

The time histories in Figs. I6l8l herein. relative to dissipative blocks and a dissipative soft layer, 
again exhibit amplification, beating phenomena, long codas and spatial variability of response in 
all three cities. The spectra exhibit so-called 'splittings' pHI, which are responsible for beating 
phenomena and part of the duration. The durations (time during the which amplitude is more 
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than 0.1 of the maximal ampHtude) are not notably different for the three cities (i.e., maximal 
duration of around 120s) for the three configurations. The entries in Table IS indicate that the 
vulnerability indices of the blocks are significantly-large even in the case of dissipative media. 
The spectra are more irregular (i.e., the number of peaks is larger) when the average separation 
of the blocks is smaller. This may constitute an indication of greater structure-soil-structure 
and/or greater block-to-block interaction responsible for more efficient excitation of quasi Love 
modes (these are neither the rigid base block modes, nor the Love modes in the absence of the 
buildings, but a combination of the two |11|V 

6 Conclusion 

The traditional 0, j2] (and even quite recent P], (S], 0) choice, concerning the analysis of 
earthquakes in cities built on soft soil as Nice and Mexico City, has been to introduce complexity 
into the substratum while leaving the free surface flat or with large-scale topographic features 
0. The buildings are not included in such models and their response is treated separately, if 
at all, using the flat ground motion as the input. This has led to response predictions that are 
more in line with what has actually been observed during tremors in a variety of cities (Nice, 
Los Angeles, Mexico City, etc.). Nevertheless, one or several of the features, namely duration, 
peak velocities and spatial variability, of observed response, differ from those of the predictions. 
This is possibly due to the fact that the flne structure of the substratum is usually unknown. 
Another plausible hypothesis is that small-scale irregularities on the free surface and on interface 
between the foundations and the soft soil, introduced by the existence of buildings in a city, may 
contribute significantly to the overall motion of the site. 

This hypothesis was shown to be tenable for: i) a non-dissipative basement underlying a 
city of ten, non-equispaced, non-equally sized, homogenized, non-dissipative blocks in [Q, and 
ii) a dissipative basement underlying a city with an infinite number of equispaced, identical, 
dissipative blocks in ^H]- Herein, we again considered cities with ten non-equispaced non-equally 
sized, homogenized blocks, but the latter, as well as the underlying layer, were considered to be 
dissipative. 

We found that the introduction of dissipation does not eliminate the previously-found anoma- 
lous effects arising from the presence of the blocks, which take the form of peak ground and block 
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motion amplifications, beatings, and long codas. As expected, the durations turned out to be 
smaller than in the non-dissipative case, but the vulnerability indices remained significantly 
large. When the spacing between blocks was reduced, we found some evidence of increased 
structure-soil-structure and block-to-block interactions, but, contrary to expectation, this did 
not have significant effects on either the duration of motion or on the peak amplification. 

These anomalous effects, as well as their spatial variability, are qualitatively the same as 
those observed in Mexico City and underline the fact that the individual buildings or groups of 
buildings (blocks) play a very active role in the overall motion of a city submitted to a seismic 
disturbance. 

In the future, it will be necessary to examine to what extent the anomalous response is 
affected by the location and type of seismic source, as well as by the duration of the pulse 
radiated from this source (note that studies such as [3], ^Hl take no account of the buildings of 
the city). 
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Figure 2: Time history of the displacement (left panel) and displacement spectrum (right panel). 
Solid-line: numerical solution, dashed-line: semi-analytical solution 
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Figure 3: Snapshot of the total (i.e., incident plus scattered) displacement field, at t = 15s, for 
flat ground underlain by a soft layer and hard half space. 




Figure 4: Snapshot of the total (i.e., incident plus scattered) displacement field, at t = 15s, for 
city with 10 blocks. 
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Figure 5: Time record of the velocity on top of building no. 6, for city with 10 blocks. Solid 
line: dissipative media, dashed line: non-dissipative media. 
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m-2=l.3416 









Figure 6: Time records of total particle velocity for 'city' with ten blocks having different 
spacing djj+i {j = 1,2, ...,9). Each row of the figure depicts the particle velocity (in s): at the 
center of the top of the j-th block (left), the center of the ground segment between the j-th 
and the j + 1-th block (middle) and at the top of the j + 1-th block (right). The solid curves 
in all the subfigures represent the particle velocity at ground level in the absence of the blocks. 
The vulnerability indices Rj at the top of the j-th block and Rjj+i on the ground between the 
j-th and the j + 1-th block, are indicated at the top of each subfigure. The abscissas designate 
time, and range from to 200s. Note that the scales of the ordinates do not vary from one 
subfigure to another. Below each time record appears the modulus of the velocity spectrum of 
the previous quantity. The abscissas designate frequency, and range from to 2Hz. 
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Figure 7: Time records of total particle velocity for 'city' with ten blocks having different 
spacing djj+i {j = 1,2, ...,9). Each row of the figure depicts the particle velocity (in s): at the 
center of the top of the j-th block (left), the center of the ground segment between the j-th 
and the j + 1-th block (middle) and at the top of the j + 1-th block (right). The solid curves 
in all the subfigures represent the particle velocity at ground level in the absence of the blocks. 
The vulnerability indices Rj at the top of the j-th block and Rjj+i on the ground between the 
j-th and the j + 1-th block, are indicated at the top of each subfigure. The abscissas designate 
time, and range from to 200s. Note that the scales of the ordinates do not vary from one 
subfigure to another. Below each time record appears the modulus of the velocity spectrum of 
the previous quantity. The abscissas designate frequency, and range from to 2Hz. 
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Figure 8: Time records of total particle velocity for 'city' with ten blocks having different 
spacing d^j+i {j = 1,2, ...,9). Each row of the figure depicts the particle velocity (in s): at the 
center of the top of the j-th block (left), the center of the ground segment between the _7-th and 
the j + 1-th block (middle) and the j + 1-th block (right). The solid curves in all the subfigures 
represent the particle velocity at ground level in absence of blocks. The vulnerability indices 
Rj at the top of the j-th block and Rj,j+i on the ground between the j-th and the j + 1-th 
block, are indicated at the top of each subfigure. The abscissas designate time, and range from 
to 200s. Note that the scales of the ordinates do not vary from one subfigure to another. 
Below each time record appears the velocity spectrum of the previous quantity. The abscissas 
designate frequency, and range from to 2Hz. 
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Bi 


B2 


B^ 


Ba 


B^ 


Bq 


B7 


Bf, 


B^ 


-Bio 


height hi 


50 


50 


50 


60 


60 


60 


70 


70 


70 


70 


width Wi 


30 


40 


50 


30 


40 


60 


30 


40 


50 


60 


spacing dj^^^ 


70 


90 


60 


80 


100 


60 


90 


80 


60 






30 


40 


30 


40 


50 


30 


40 


40 


30 




"-1,1+1 


10 


20 


10 


20 


30 


10 


20 


20 


10 





Table 1: Geometrical parameters of the ten-block idealized cities. The units are meters 







R, 










R- 






-Rio 




9.2 


10.2 


7.5 


12.6 


9.4 


6.1 


10.9 


7.9 


6.9 


6.5 




7.5 


7.1 


4.0 


6.7 


5.9 


3.2 


6.0 


5.4 


4.0 


3.3 




4.2 


9.5 


5.9 


5.6 


7.6 


4.0 


5.0 


6.6 


4.0 


2.7 


C3 


2.6 


6.2 


4.9 


3.8 


6.1 


4.8 


5.7 


7.5 


4.6 


2.7 



Table 2: Vulnerability indices, (bold for non-dissipative media) 
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